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A study of event shape observables in proton-antiproton collisions at y/s=1.96 TeV is presented. 
The data for this analysis were recorded by the CDF II detector at the Tevatron collider. The 
variables studied are the transverse thrust and thrust minor, both defined in the plane perpendicular 
to the beam direction. The observables are measured using energies from unclustered calorimeter 
cells. In addition to studies of the differential distributions, we present the dependence of event 
shape mean values on the leading jet transverse energy. Data are compared with pythia Tune 
A and to resummed parton level predictions that were matched to fixed order results at NLO 
accuracy (NLO+NLL). Predictions from pythia Tune A agree fairly well with the data. However, 
the underlying event contributes significantly to these observables, making it difficult to make direct 
comparisons to the NLO+NLL predictions, which do not account for the underlying event. To 
overcome this difficulty, we introduce a new observable, a weighted difference of the mean values 
of the thrust and thrust minor, which is less sensitive to the underlying event, allowing for a 
comparison with NLO+NLL. Both pythia Tune A and the NLO+NLL calculations agree well within 
the 20% theoretical uncertainty with the data for this observable, indicating that perturbative QCD 
successfully describes shapes of the hadronic final states. 



I. INTRODUCTION 



* Deceased 

tWith visitors from "University of Massachusetts Amherst, 
Amherst, Massachusetts 01003, ^Istituto Nazionale di Fisica Nu- 
cleare, Sezione di Cagliari, 09042 Monscrrato (Cagliari), Italy, 
c University of California Irvine, Irvine, CA 92697, ''University of 
California Santa Barbara, Santa Barbara, CA 93106 e University 
of California Santa Cruz, Santa Cruz, CA 95064, f CERN,CH- 
1211 Geneva, Switzerland, ^Cornell University, Ithaca, NY 14853, 
''University of Cyprus, Nicosia CY-1678, Cyprus, 'University Col- 
lege Dublin, Dublin 4, Ireland, 3 University of Fukui, Fukui City, 
Fukui Prefecture, Japan 910-0017, k Univcrsidad Ibcroamcricana, 
Mexico D.F., Mexico, 'Iowa State University, Ames, IA 50011, 
m University of Iowa, Iowa City, IA 52242, "Kinki University, 
Higashi-Osaka City, Japan 577-8502, "Kansas State University, 
Manhattan, KS 66506, p Univcrsity of Manchester, Manchester M13 



The hadronic final states produced in hard collisions 
can be characterized by a number of variables that de- 



9PL, England, 9 Queen Mary, University of London, London, El 
4NS, England, r Muons, Inc., Batavia, IL 60510, "Nagasaki In- 
stitute of Applied Science, Nagasaki, Japan, *National Research 
Nuclear University, Moscow, Russia, "University of Notre Dame, 
Notre Dame, IN 46556, " Universidad de Oviedo, E-33007 Oviedo, 
Spain, ""Texas Tech University, Lubbock, TX 79609, x Universidad 
Tecnica Federico Santa Maria, llOv Valparaiso, Chile, "Yarmouk 
University, Irbid 211-63, Jordan, "On leave from J. Stefan Insti- 
tute, Ljubljana, Slovenia, 



4 



scribe the distribution of outgoing particles in the event. 
These are referred to as event-shape variables. Mea- 
surements of these variables in e + e~ and deep-inelastic 
scattering experiments [1] allowed a determination of 
the strong coupling constant a s and its renormalization 
group running Q , color factor fits of the quantum chro- 
modynamics (QCD) gauge group !3j, and more recently, 
studies of non-perturbative corrections to QCD reviewed 
in [H . The goal of the analysis here is to measure event 
shapes in proton-antiproton collisions and to study their 
dependence on the transverse energy of the leading jet. 
The data are compared with pythia Tune A [j| Q and 
to resummed next-to-leading-logarithm parton level pre- 
dictions that were matched to fixed order results at next- 
to-leading-order accuracy (referred to as NLO+NLL). 
This study contributes to our understanding of the un- 
derlying event (UE) in a hard-scattering process, and 
illustrates the need to include underlying event contri- 
butions when comparing data with perturbative QCD in 
hadron-hadron collisions. 

In general, event-shape variables describe geometric 
properties of the energy flow in QCD final states. They 
are related to jet-finding algorithms, which are used to 
categorize events according to their topology. However, 
they differ significantly in that event-shape variables en- 
code information about the energy flow in the overall 
event. A single parameter can describe, for example, 
the transition between a configuration with all particles 
flowing along a single axis and a configuration where the 
energy is distributed uniformly over solid angle. Event- 
shape variables therefore provide an alternative way to 
characterize an event compared to others based on jet- 
finding algorithms. Furthermore, they have the advan- 
tage of being free of the arbitrariness associated with jet 
definition, i.e., being either cone or cluster in type, cone 
sizes, splitting/merging fractions, etc. 

The earliest studies of event shapes in hadron-hadron 
collisions were performed at the ISR |8| and the SPS 

in the late 1970s and focused on tracing the emer- 
gence of jet-like structures. A decade later event shapes 
were measured at a hadron-hadron collider during Run 

1 at the Tevatron, where variants of jet broadening and 
thrust were measured by CDF [Io| and DO [HI- Ab- 
sent from all of these studies was a direct comparison 
to perturbative QCD. This was in large part due to the 
intrinsic theoretical difficulties associated with modeling 
of the hadron collider environment. However, recently 
a number of tools for investigating event-shape variables 
beyond leading order have been developed [7(], allowing 
for comparisons with hadron collider data. 

In this paper, we report a measurement of event-shape 
variables using energies from unclustered calorimeter 
cells in events with transverse energy of the highest- 
energy (leading) jet ranging from 100 to 300 GeV. 
Events were produced at the Tevatron collider in proton- 
antiproton collisions at a center-of-mass energy of 1.96 
TeV and were recorded by the CDF II [l2j detector. The 
data used in this analysis were collected from February 



2002 to August 2004, with an integrated luminosity of 
385 pb" 1 . Section II contains a brief description of the 
NLO+NLL theoretical predictions used in this analysis. 
In Sec. Ill the observables of interest are introduced 
and the effects of hadronization and the underlying 
event are examined. Section IV outlines the features of 
the CDF II detector relevant for this analysis. Event 
pre-selection and the reconstructed objects used in the 
analysis are described in Sec. V. Section VI reviews 
the measurement, describes instrumental effects on the 
measurements, and explains how the measurements 
are corrected for these effects. Section VII lists the 
sources of systematic uncertainties that affect the final 
results. Finally, results are presented and summarized 
in Sees. VIII and IX, respectively. 



II. NLO+NLL THEORY 

The calculation of event-shape-variable distributions 
at the parton level in perturbative QCD is divided into 
two regimes: fixed order and resummed calculations [Zj. 
Event shape observables considered in this analysis have 
the property that large values of the variable coincide 
with the emission of one or more hard partons at large 
angles relative to the parent parton. In this regime, the 
distribution is well described by a traditional perturba- 
tive expansion in powers of the strong coupling, a s . This 
method provides an accurate description over most of the 
range of the variables, e.g. transverse thrust discussed 
later in this paper. 

The method breaks down for small values of the event- 
shape variables. In fact, all fixed order calculations di- 
verge in the limit that the event shape variable goes to 
zero. In this region, the differential cross section is pri- 
marily sensitive to gluon emissions that are soft compared 
to the hard scale of the event and/or collinear with one 
of the hard partons. Such radiation has relatively large 
emission probabilities due to logarithmic enhancements. 
In this case, each power of a s in the perturbative expan- 
sion is accompanied by a coefficient that grows as In 2 X/y, 
where y is the variable of interest. This enhances the im- 
portance of higher order terms in the series and the naive 
requirement that a s be small is no longer sufficient to ren- 
der these terms negligible. To obtain meaningful predic- 
tions in the region y — ^ it is necessary to perform an all- 
orders resummation of the enhanced logarithmic terms, 
which can be performed with next-to-leading-logarithm 
(NLL) precision. 

The parton level theoretical predictions used in this 
analysis correspond to fixed-order results at next-to- 
leading-order accuracy matched to resummed results at 
next-to-leading-logarithm accuracy. Theoretical fixed or- 
der results are obtained from the Monte-Carlo integra- 
tion program nlojet-) — h [l3j], while the resummed re- 
sults arise from the Computer Automated Expert Semi- 
Analytical Resummer (CAESAR) 0, These the- 
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oretical calculations include initial and final-state radia- 
tion, but do not include multiple parton interactions or 
beam remnant models. They all use CTEQ 6.1 [l6| par- 
ton distribution functions (PDFs). 

Another technical restriction of NLL calculations is 
that they apply only to "global" variables (i.e., are sen- 
sitive to radiation in all directions). This requirement is 
in direct conflict with the realities of any collider exper- 
iment; namely the limited detector coverage at large ra- 
pidities. However, the variables considered in this paper 
are defined exclusively in the transverse plane perpen- 
dicular to the beam axis. Therefore, for sufficiently large 
values of the maximum accessible rapidity, the contribu- 
tion from the excluded kinematic region is expected to 
be small and the full global predictions for the studied 
variables should remain valid for ln(y) < krj ma x-, where 
Vmax is the maximum detector pseudorapidity coverage 
and k is a constant dependent on the variable y Q. 

III. DEFINITIONS OF THE VARIABLES 

A. Transverse Thrust and Thrust Minor 

The transverse-thrust variable, in analogy to the usual 
thrust, is defined as [l(: 

71 

T± = max ^-jj , (1) 

£l?r,*l 
»=i 

where q± + is the transverse momentum of the i'th ob- 
ject, where the object is either an outgoing parton, an 
outgoing stable particle, or a calorimeter cell. The sum 
runs over all objects in the final state and the thrust axis 
ut is defined as the unit vector in the plane perpendic- 
ular to the beam direction that maximizes this expres- 
sion. For an event with only two back-to-back outgoing 
objects T± = 1. In the case of a perfectly cylindrically 
symmetric event the transverse thrust takes on the value 
Tj_ = 2/ir w 0.637. Historically the majority of event- 
shape variables are defined so that they vanish in the 
limit of two back-to-back objects. Therefore it is conve- 
nient to define r = 1— Tj_, which has this property. Here- 
after, any discussion of the variable called thrust shall 
refer to the quantity r, where < r < 1 — 0.637. 

The thrust axis fix and the beam direction z together 
define the event plane in which the primary hard scat- 
tering occurs. The transverse thrust minor is defined as: 

n 

^ \q±,i ■ n m \ 

1=1 



where n m = ut x z. The observable T m i n is a measure of 
the out-of-plane transverse momentum and varies from 
zero for an event entirely in the event plane to 2/ir for a 
cylindrically symmetric event. 

It should be noted that the authors of Ref. 0] also pro- 
posed an alternative definition for event-shape variables 
at hadron colliders to specifically deal with the issue of 
limited detector coverage. As originally envisioned, the 
event-shape variables were to be defined over outgoing 
objects in a reduced central region and rendered "indi- 
rectly" global by the addition of a "recoil" term event- 
by-event. The "recoil" term is defined by momenta in 
the same central region, but would introduce an indirect 
sensitivity to momenta outside the region. The proposed 
recoil term is essentially the vector sum of the transverse 
momenta in this central region (which by conservation of 
momentum is equal to the vector sum of transverse mo- 
menta outside the region). However, Monte Carlo studies 
done for this analysis showed that there was almost no 
correlation between the event-shape variables and the re- 
coil term. As a result, this alternative definition was not 
pursued. 

A number of other event-shape variables (broadenings, 
hemisphere masses, etc.) whose definitions include lon- 
gitudinal components of the final state momenta have 
also been proposed Q. Our studies indicate that these 
variables are very sensitive to detector resolution in the 
forward region (see Sec. IV). Hence, we focus on a study 
of the observables r and T m i n defined in the transverse 
plane. By construction these two quantities are infrared 
and collinear safe. 



B. Hadronization and the Underlying Event 

The NLO+NLL parton level calculations do not in- 
clude hadronization effects and they do not have a model 
of the underlying event. In particular, they do not in- 
clude beam-beam remnants nor multiple parton interac- 
tions [13. We use pythia 6.1 5] with CTEQ5L PDF 
sets [18| to study the effects of hadronization and the 
underlying event on the transverse thrust, r, and thrust 
minor, T m i n . The underlying event corresponds to par- 
ticles that arise from the beam-beam remnants or from 
multiple-parton interactions. Figure [1] shows a compar- 
ison of the distributions of transverse thrust and thrust 
minor between the NLO+NLL parton level theoretical 
predictions, and PYTHIA without underlying event at the 
parton level and after hadronization. The comparison 
is made for events with leading jet transverse energy 

E lead.jet > 2QQ ^y. ^ trangverse 

energy is defined in 
Sec. IV. The plot shows that pythia and the NLO+NLL 
parton level predictions have similar shapes for both the 
transverse thrust and the thrust minor. However, for the 
transverse thrust, the PYTHIA distribution is shifted to- 
ward larger values over the entire range of the variable. 
Furthermore, hadronization in pythia produces only a 
small shift of the event-shape distributions towards val- 
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ues larger than pythia without hadronization, a result 
expected from LEP [Hj]). 



Leading Jet E T > 200 GeV 




FIG. 1: Predictions of the transverse thrust and thrust minor 
distributions for Rf 1 J e< greater than 200 GeV from a parton 
level NLO+NLL calculation and from pythia without an un- 
derlying event at the parton level and without an underlying 
event at the hadron level (i.e. after hadronization). 

Figure [2] shows a comparison of the event-shape distri- 
butions between the NLO+NLL parton level predictions, 
pythia without underlying event, and pythia Tunc A. 
pythia Tune A includes a model of the underlying event, 
which was tuned to fit the CDF Run I underlying event 
data. We see that the underlying event not only shifts 
the means towards higher values, but also significantly 
affects the overall shape of the distributions. Figure [3] 
shows mean values as a function of leading jet trans- 
verse energy. There is very little difference in the mean 
values for pythia Tune A at the parton and hadron lev- 
els. The additional partons from multiple-parton inter- 



actions saturate the event-shape variable distributions to 
a point where the "re-shuffling" of momenta that occurs 
at hadronization has little effect on the variable. 




FIG. 2: Predictions of the transverse thrust and thrust mi- 
nor distributions for E^ a ' 3ei greater than 200 GeV from a 
parton-level NLO+NLL calculation and from pythia at the 
parton level without an underlying event and at the parton 
level with an underlying event (Tune A). 

We conclude that the underlying event significantly af- 
fects the distributions of t and T m i n . As a result, a di- 
rect comparison of event-shape variable distributions in 
data with the NLO+NLL parton level predictions is not 
possible. However, a quantity less dependent on the un- 
derlying event can be constructed from the average val- 
ues of the thrust and thrust minor. The dependence of 
this quantity on the leading jet transverse energy might 
then allow for a more meaningful comparison between 
NLO+NLL parton level predictions and the measured 
data. To this end we begin by considering the definitions 
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of the thrust Eq. (JXJ) and thrust minor Eq. ([2]). Sepa- 
rating the final state into hard (q l ^ ard ) and soft (g^*) 
components and recognizing that the thrust axis is de- 



termined almost entirely by the hard component, we see 
that the transverse thrust and thrust minor can be writ- 
ten approximately as: 



^T^-max^ricos^l y q ^f*(i _ | cos ^/*| 



rir — / j 



(3) 



Tn 



E <?i ard | sin (f> hard | E qS ° ft I sin ^° ft I 



(4) 



where (p hard and SO ^* represent the angle between the 
thrust axis and the hard and soft components, respec- 
tively. The soft underlying event is expected to be 
on average uniform over the transverse plane, therefore 



1 — t so ^ 1 sa T^/J f=3 2/tt. An expression whose numer- 
ator is less dependent on the underlying event can be 
constructed by taking a weighted difference between the 
mean values of the thrust and thrust minor as follows: 



- hard \cos<b hard \ 



I y q 'r d I sin j> hard I \ / E & ard - m - ax E 



where a — 1 — 2/tt, /? = 2/7r, and the sums run over 
objects from the hard scattering ("hard") and from the 
soft underlying event ( "soft" ) . The soft underlying event 
in this expression is in the denominator where its con- 
tribution is overshadowed by the hard-scattering term. 
Furthermore, an additional correction factor, Jmc-> can 
be computed from pythia tune A 0, Q generated with 
and without multiple parton interactions: 



IMC 



EK oMPJ l+EI^ 

El NoMPI I 
\9± 



WithMPI I 



Now we define a new variable: 

D((t) , (T mm )) = jMc{a (T min ) - /3 (r)), 



(6) 



(7) 



which should be less dependent on the underlying event. 
The 7mc correction is applied to the data for compar- 
isons to theory. We call D((t) , (T m i n )) the "thrust differ- 
ential". It should be insensitive to the underlying event 
activity in the event, thereby allowing more meaningful 
comparisons to perturbative calculations of event-shape 
variables. The variable ranges between and 0.1 and 
vanishes in both limiting cases of cylindrically symmet- 
ric and pencil-like events. The variable allows probing 
the relative contributions of pQCD and non-pQCD pro- 
cesses to the distributions of the event-shape variables. 



It is important to note that the Jmc correction factor 
differs from unity by no more than 13% over the range 
of the leading jet E^ a ' Je * threshold, as shown in Fig. |U 
Also, it can be seen in Fig. [5] that MC simulations show 
a strongly-reduced effect of the underlying event on the 
thrust differential relative to its effect on the transverse 
thrust and thrust minor shown in Fig. [3] Experimen- 
tal measurements of transverse thrust and thrust minor, 
of the thrust differential, and of the dependence of the 
latter quantity on the transverse energy threshold of the 
leading jet are presented in the next sections. 



IV. CDF II DETECTOR 

Data used in this analysis were recorded with CDF II, 
a general-purpose detector, designed for precision mea- 
surements of the energy, momentum, and trajectories of 
particles produced in proton-antiproton collisions. This 
section provides a brief overview of the components rele- 
vant to our analysis. A detailed description of the entire 
detector can be found elsewhere (T^ |. 

CDF II uses a spherical coordinate system with the z 
axis oriented along the proton beam direction and az- 
imuthal angle <f> measured around the beam axis. The 
polar angle 6 is measured with respect to the positive 
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FIG. 3: Predictions of the mean values of the transverse 
thrust and thrust minor as a function of the leading jet trans- 
verse energy from a parton level NLO+NLL calculation and 
from pythia at the parton level without an UE (pythia Par- 
ton), at the parton level with an underlying event (Tune A 
Parton), and at the hadron level with an underlying event 
(Tune A Hadron). 
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FIG. 4: Mean value of 7mc, from Eq. (6), as obtained from 
pythia at the parton level with and without multiple parton 
interactions. 
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FIG. 5: Predictions of the mean values of the thrust differen- 
tial as a function of the leading jet transverse energy from a 
parton level NLO+NLL calculation and from pythia at the 
parton level with (Tune A) and without an underlying event 
(Parton). The plot indicates that the underlying event has 
only a small effect on the mean value of the thrust differen- 
tial. 



z (proton-beam) direction used to define pseudorapidity 
r/ = - In [tan(f)]. 

The CDF II tracking system is used for reconstruc- 
tion of primary interaction vertices and particle tracks, 
and is placed inside a 1.4 T solenoidal magnet. An in- 
ner, single-sided silicon microstrip detector (Layer 00) is 
mounted directly on the beam pipe, at an inner radius of 
1.15 cm and an outer radius of 2.1 cm. A five-layer sili- 
con microstrip detector (SVX II) is situated at the radial 
distance of 2.5 to 11 cm from the beam line, and con- 
sists of three separate barrel modules with a combined 
length of 96 cm. Three of the five layers combine a r- 
<f> measurement with a z-coordinate measurement while 
the remaining two layers combine r-<f> with small- angle 
stereo views 1.2°. Three additional intermediate silicon 
layers (ISL) are positioned between 19 and 30 cm from 
the beam line. The silicon tracker is surrounded by the 
central outer tracker (COT), an open-cell drift chamber 
providing up to 96 measurements of a charged particle 
track over the radial region from 40 to 137 cm. The 
pseudorapidity region covered by the COT is \rj\ < 1.0. 

The CDF II tracking system is surrounded by elec- 
tromagnetic and hadronic calorimeters, whose cells are 
arranged in a projective tower geometry. The central 
electromagnetic (CEM), central hadronic (CHA), and 
wall hadronic calorimeters consist of lead (electromag- 
netic) and iron (hadronic) layers interspersed with scin- 
tillator. The pseudorapidity region covered by these 
calorimeters is \rj\ < 1.3. The segmentation of the cen- 
tral calorimeters is 15° in (f> and 0.1 units in r\ and again 
15° in (f> but 0.2 to 0.6 in 77 in a forward "plug" region. 
The measured energy resolutions for the CEM and CHA 
are a{E T )/E T = 13.5%/ v / ^©2% and a(E T )/E T = 
75%/v^r © 3%, respectively. Here E T = Esin9 is the 
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transverse energy deposited in a calorimeter tower. En- 
ergies are measured in GeV. Additional calorimetry ex- 
tends the coverage in the forward direction to \rj\ < 3.6. 
The forward electromagnetic calorimeter is constructed 
of lead and scintillator layers with an energy resolution 
of a{E T )/E T = 16%/vW © 1%. The forward hadronic 
calorimeter is made of iron and scintillator layers with an 
energy resolution of a(E T )/E T = 80%/y/E^ © 5%. 

V. EVENT SELECTION 

A. Triggers 

Events were collected using single-jet triggers with Et 
thresholds of 50 (J050), 70 (J070), and 100 (J100) GeV. 
Pre-scale factors are applied to J050 and J070 jet triggers 
so as not to saturate the available trigger bandwidth; 
typical values of pre-scale factors are 8 and 50 for J050 
and J070, respectively. The J100 trigger is not prescaled. 

B. Jet reconstruction algorithm 

While event shape variables are calculated from un- 
clustered calorimeter cell energies, analyzing their depen- 
dence on leading-jet Et requires use of a jet-energy algo- 
rithm for the leading jet scale. Based on calorimeter in- 
formation, jets are reconstructed using a cone algorithm 
[20| . The algorithm starts with the highest transverse 
energy tower and forms preclusters from an unbroken 
chain of continuous seed towers with transverse energies 
above 1 GeV within a window of 7 x 7 towers centered 
on the originating seed tower. If a seed tower is out- 
side this window, it is used to form a new precluster. 
The coordinates of each precluster are the .Ei-weighted 
sums of <p and r\ of the seed towers within this precluster. 
In the next step, all t owers with Et > 0.1 GeV within 
R = ^{A(f>) 2 + (At]) 2 = 1.0 of the precluster are merged 
into a cluster, and its {rj, 0)-coordinates are recalculated. 
This procedure of calculating cluster coordinates is iter- 
ated until a stable set of clusters is obtained. A cluster 
is stable when the tower list is unchanged from one itera- 
tion to the next. If the clusters have some finite overlap, 
then an overlap fraction is computed as the sum of the 
transverse energies of the common towers divided by the 
Et of the smaller cluster. If the fraction is above a cutoff 
value of 0.75, the two clusters are combined. If the frac- 
tion is less than the cutoff, the shared towers are assigned 
to the closer cluster. The raw energy of a jet is the sum of 
the energies of the towers belonging to the correspond- 
ing cluster. Corrections are applied to the raw energy 
to compensate for the non-linearity and non-uniformity 
of the response of the calorimeter, the energy deposited 
inside the jet cone from sources other than the assumed 
leading jet-parent parton, and the leading parton energy 
deposited outside the jet cone. A detailed description of 
this correction procedure can be found in 21]. 



TABLE I: Summary of the data samples, trigger paths, and 
number of events present after the offline selection criteria. 



E le.ad. je t 


Trigger 


Number of Events 


100 


J050 


52546 


150 


J070 


17850 


200 


J100 


26207 


300 


J100 


3126 



C. Offline selection 

Cosmic ray events are rejected by applying a cutoff 
on the significance of the missing transverse energy Jf*T 
[H], defined as $ T where T,E T = T, t E T is the 
total transverse energy of the event, as measured using 
calorimeter towers with E T above 100 MeV. The cutoff 
values are 5.0, 6.0, and 7.0 GeV 1 / 2 for data collected 
using jet triggers with thresholds of 50, 70, and 100 
GeV, respectively. 

To ensure fully efficient vertex and track reconstruc- 
tion, we require events with a single interaction as evi- 
denced by having only one reconstructed primary inter- 
action vertex with | z \ < 60 cm. 

Only events with both leading jets in the central region 
(|?7| < 0.7) are selected. Events are categorized accord- 
ing to the transverse energy threshold of the leading jet 
as in Ref (tJ- The event categories, trigger paths, and 
number of events after selection criteria are summarized 
in Table H 



VI. MEASUREMENT AND INSTRUMENTAL 
UNCERTAINTIES 

The measurement of event-shape variables is per- 
formed with unclustered calorimeter towers over the de- 
tector's full rapidity range < 3.5. Towers used in the 
measurement are required to have minimum E^ of 100 
MeV. 

In general, measurement of the event-shape variables 
will be distorted by instrumental effects and a correction 
factor is needed to account for this. Figure [5] shows the 
dependence of D((t), (TM in )) on the leading jet trans- 
verse energy for pythia Tune A at the hadron level and 
at the calorimeter level after full CDF detector simula- 
tion. While the detector effects do not appreciably affect 
distributions of r and T TO ,- n , they induce a significant sys- 
tematic shift of the thrust differential. Sources for this 
shift have been studied, and the results are given below. 

As a result of the magnetic field, the energy flow of 
an event as measured by the calorimeter will be broader 
than in the absence of the field. To estimate the magni- 
tude of this effect on the thrust differential, MC simulated 
particles at the hadron level were propagated to the first 



10 



ter granularity. 



qO.04 
0.038 
0.036 
0.034 
0.032 

0.03 
0.028 
0.026 
0.024 
0.022 

0.02 



Tune A Hadron 
Tune A + CDFSim 



100 



150 



200 



250 300 

E tead.]et {GeV j 



FIG. 6: The effect of CDF detector simulation on the thrust 
differential as a function of the leading jet transverse energy. 



active layer of the calorimeter under the influence of a 
1.4 T B-field. The direction of each particle was calcu- 
lated from the z coordinate of the primary interaction 
vertex and the point of impact on the first calorimeter 
layer. The effect is found to be ~ 2% of the values of the 
thrust differential and negligible compared to the effect 
of calorimeter granularity described below. 

To estimate the effect on the thrust differential of the 
calorimeter energy resolution, we smear the energy of the 
particles in the MC simulation according to a gaussian 
with the lcr resolution quoted in section IV. Smearing 
changes the thrust differential by < 1% of its value. 

Turning now to the effect of the calorimeter granu- 
larity, we note that when a particle above threshold is 
detected, the location returned by the system is the cen- 
ter of the tower and not the exact location of the shower 
within the tower. As a result, there is an error associ- 
ated with the granularity of the calorimeter. In order to 
understand this effect on the thrust differential, the seg- 
mentation of the calorimeter is imposed on MC simulated 
particles at the hadron level. 

The effects of these instrumental errors on the thrust 
differential are shown in Fig. [7] The granularity of the 
calorimeter is the primary source of the shift of this vari- 
able. The shift is taken into account by a bin-by-bin 
correction to the thrust and thrust minor distributions, 
which is propagated to D((t), (TWin))- 

In the model calculations referenced here, event-shape 
variables are defined over all particles in the final state, 
including those with arbitrarily small momenta. In order 
to understand how a cut on the transverse energy affects 
the variables, we vary the threshold on towers from 
100 MeV (default) through 200 and 300 MeV. Figure M 
shows that the thrust differential is rather insensitive to 
the cut on transverse energy at low q±. While the dis- 
tributions of thrust and thrust minor get narrower with 
increasing j^ ad -i et threshold, the effect on the thrust dif- 
ferential is negligible compared to the effect of calorime- 
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FIG. 7: Simulation of individual instrumental effects on the 
thrust differential as a function of the jet energy threshold. 
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FIG. 8: Effect of tower Et threshold on the thrust differential 
as a function of the leading jet transverse energy threshold. 



VII. SYSTEMATIC UNCERTAINTIES 

The sensitivity of the thrust differential to various un- 
certainties in the event selection procedure is evaluated as 
follows. For each source of systematic uncertainty, a "de- 
fault" and "deviated" variable is constructed. The "de- 
fault" variable is the result of the standard set of cuts de- 
fined earlier in this paper, while the "deviated" variable 
is the result of varying a particular parameter within its 
uncertainty. For each of the data samples corresponding 
to the four different values of transverse energy thresh- 
olds of the leading jet, the systematic uncertainties on 
the thrust differential are calculated as the difference be- 
tween the "deviated" and the "default" results. Each 
individual source of uncertainty is then added in quadra- 
ture to the statistical uncertainty of each data point. 
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A. Jet Energy Scale 

Theoretical predictions for the event-shape variables 
are parametrized as functions of the leading jet transverse 
energy. To evaluate the uncertainty on the leading jet 
transverse energy due to the jet energy corrections, we 
use a parametrization that under- and over-estimates the 
leading jet energy by one standard deviation in the jet 
energy scale [2l[ and then re-run our event selection. The 
difference between the default and the deviated variable 
is assigned as a systematic uncertainty. As expected, 
most of the jet energy scale error cancels in the ratio of 
sums of Eqs. (1) and (2) and in the calculation of D, 
hence the resulting systematic uncertainty is small. 

B. Detector Hermeticity 

The primary interaction vertex is required to lie within 
60 cm from the center of the detector in order to en- 
sure that the majority of the event is contained within 
the detector. This analysis of event-shape variables uses 
calorimeter information in the far forward regions of the 
detector. As a result, the further a collision occurs from 
the nominal interaction point the greater the possiblity 
that particles fall beyond the detector's coverage. To 
evaluate the uncertainty due to this effect we require 
a tighter cut on the z position of the primary vertex, 
\z\ < 20 cm. The difference in the variable between the 
default and the tight cut is then assigned as a systematic 
uncertainty. 

C. Effect of Pile-up 

In the event selection we specifically require events 
with a single vertex; however, it is possible that two ver- 
tices that lie very close to each other are reconstructed as 
a single vertex. This "pile-up" effect is especially likely 
at high values of the instantaneous luminosity. To evalu- 
ate the uncertainty due to this effect we separate events 
in each data sample into high (average of 3 primary in- 
teractions per bunch crossing) and low (average of 1.5 
primary interactions per bunch crossing) luminosity sub- 
sets with approximately equal numbers of events. The 
thrust differential is then compared between subsets and 
the difference is taken as a measure of the systematic 
uncertainty. 

A summary of all uncertainties affecting this measure- 
ment are given in Table HH together with the values of 
the thrust differential. 



VIII. RESULTS 

The distributions of the transverse thrust and thrust 
minor, uncorrected for detector effects, are shown in 



Fig.Ufor events with the E l ^ addet greater than 200 GeV. 
Distributions for other fi T ™ ' Jef thresholds can be found 
in [23[. There is not much difference between pythia 
Tune A at the hadron level and the detector level (CDF- 
SIM) . Tune A describes the data fairly well although not 
perfectly. The distribution of the thrust minor is slightly 
broader than the Tune A prediction (i.e. there is slightly 
more energy out of the plane than predicted by Tune A) . 
The parton level NLO+NLL predictions deviate signifi- 
cantly from the data since they have no underlying event. 
For events with leading jet transverse energy greater than 
200 GeV, the mean value of the r distribution shifts from 
0.039 ± 0.001 to 0.070 ± 0.001 (parton level NLO+NLL 
to experiment), while the RMS remains unchanged at 
0.040 ± 0.001. The mean value of the T m i n distribution 
shifts from 0.142 ± 0.002 to 0.206 ± 0.002 with its RMS 
decreasing from 0.099 ± 0.001 to 0.087 ± 0.001. 

Figure [TO] shows the thrust differential as a function of 
j^iead.jet j n ^ g ^he data have been corrected for de- 
tector effects. The data are compared with pythia Tune 
A and with parton level NLO+NLL calculations. The 
NLO+NLL predictions shown in this figure correspond 
to a particular choice of renormalization and factoriza- 
tion scale, namely the transverse energy of the leading 
jet; the theoretical uncertainty on thrust differential is 
approximately 20% @, HH (theoretical uncertainties on 
transverse thrust and thrust minor are smaller and are of 
the order of 10% or less [24|). By construction this ob- 
servable all but eliminates the sensitivity to the underly- 
ing event. Based on 7a/c variation studies we estimated 
that the residual effect of the UE is less than few percent. 
Both pythia Tunc A and the NLO+NLL calculations 
agree fairly well with the data, indicating that the non- 
perturbative effects are small. The corrected data and 
their associated uncertainties are listed in Table [TTJ 

IX. SUMMARY 

Event-shape variable distributions are studied using 
unclustered calorimeter energies in proton-antiproton 
collisions at a center-of-mass energy of 1.96 TeV. The 
measurements were performed using individual calorime- 
ter towers with a transverse-energy threshold of 100 MeV. 
The data are compared to pythia Tune A and to re- 
summed parton level predictions that were matched to 
fixed order results at NLO accuracy (NLO+NLL). Both 
the thrust and thrust minor distributions are sensitive to 
the modeling of the underlying event. The pythia Tune 
A distributions of the observables reproduce the experi- 
mental distributions fairly well, although not perfectly. 
The data show slightly more energy out of the hard- 
scattering plane than predicted by Tune A. These observ- 
ables can be used to improve the modeling of the underly- 
ing event. The NLO+NLL predictions differ significantly 
from both the data and from pythia Tune A since these 
calculations do not incorporate either hadronization or 
the underlying event. 
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TABLE II: Summary of the experimental values of the thrust differential and of its uncertainties. 
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Jet Energy 
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100 


315.3xl0~ 4 
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-4 


0.5xl0" 4 
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4 


1.3x10" 


-4 
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275.6 xl0~ 4 


5.9x10" 
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3.4x10" 


-4 


6.1xl0~ 4 


2.8xl0~ 4 
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-4 


4.7x10" 


-4 


7.7xl0~ 4 


7.2xl0~ 4 



Leading Jet E T > 200 GeV 

NLO+NLL 

— Tune A Hadron 
» sT Tune A + CDF Sim. 




0.04 r 



100 



-#- NLO+NLL 
■ Tune A Hadron 
T Data (Unfolded) 




150 



200 



250 300 

Ef d -i e, (GeV) 



Leading Jet E T > 200 GeV 

NLO+NLL 

■ Tune A Hadron 

Tune A + CDF Sim. 
Data 




FIG. 9: The uncorrected CDF distributions of transverse 
thrust and thrust minor for leading jet transverse energy 
greater than 200 GeV. The experimental results are compared 
with a parton level NLO+NLL calculation and with pythia 
at the hadron level (Tune A Hadron) and at the detector level 
(i.e. after CDFSIM). 



A new variable, called thrust differential, is introduced. 
It is a weighted difference of the mean values of the thrust 
and thrust minor over the event sample. By construction 
it is less sensitive to the underlying event and hadroniza- 
tion effects. Both pythia Tune A and the NLO+NLL 
calculations succeed in describing the data on the thrust 
differential. This observable allows a comparison with 
the NLO+NLL calculations, and data and theory are 
found to agree well within the 20% theoretical uncer- 



FIG. 10: The CDF corrected results for the dependence of 
the thrust differential on the transverse energy of the leading 
jet. The experimental results are compared with a parton 
level NLO+NLL calculation and with pythia Tune A at the 
hadron level. The error bars correspond to statistical and 
systematic uncertainties added in quadrature. 



tainty. 
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